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Abstract. The purely plaquette 3d Ising Hamiltonian with the spins living at the 
vertices of a cubic lattice displays several interesting features. The symmetries of the 
model lead to a macroscopic degeneracy of the low-temperature phase and prevent 
the definition of a standard magnetic order parameter. Consideration of the strongly 
anisotropic limit of the model suggests that a layered, “fuki-nuke” order still exists and 
we confirm this with multicanonical simulations. The macroscopic degeneracy of the 
low-temperature phase also changes the finite-size scaling corrections at the first-order 
transition in the model and we see this must be taken into account when analysing our 
measurements. 
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The 3d plaquette Ising Hamiltonian, with the Ising spins a 
of a 3d cubic lattice, 

H ~ ^ ^ CTiCjCf^Ci 


±1 sited at the vertices 


( 1 ) 


sits at the k = 0 limit of a one-parameter family of 3d “gonihedric” Ising Hamiltonians 
IIIIIE]- In general these contain nearest neighbour {i,j), next-to-nearest neighbour 
{{i,j)) and plaquette interactions [i,j,k,l], 

= -2K^CTiCTj + , (2) 

d.i) ((*.t>> [ipMA 

which have been hne tuned to eliminate the bare surface tension. While not gauge 
theories, the Hamiltonians are still highly symmetric. In the general case when k 7 ^ 0, 
parallel, non-intersecting planes of spins may be flipped in the ground state at zero 
energy cost, which gives a 3 x 2^^ ground-state degeneracy on an L x L x L cubic 
lattice. A low-temperature expansion reveals that this is broken at hnite temperature 
[1]. However, when k = 0, planes of spins (including intersecting ones) may be flipped 
throughout the low-temperature phase, giving a macroscopic low-temperature phase 
degeneracy of 2^^. This in turn leads to non-standard corrections to hnite-size scaling 
at the first-order transition in the model mum- The same degeneracy precludes the 
use of a standard magnetic order parameter such as m = the plaquette 

model since it would be zero throughout the low-temperature phase. 

In this paper we identify a class of suitable order parameters for the 3d 
plaquette gonihedric Ising model and using high-precision multicanonical simulations 
we investigate the scaling properties of the order parameters and their associated 
susceptibilities. 


2. Model and Observables 


If we allow for anisotropic couplings the Hamiltonian of the 3d Ising model with purely 
plaquette interactions may be written as 

L L L 


hfaniso 




x=l y=l z=l 




^x^y^z^x+l^y^z^x+l^y^z+l^x^y^z-\-l 


x=l y=l z=l 


( 3 ) 


x=l y=l z=l 


where we have indicated each site and directional sum explicitly, assuming we are on 
a cubic L X L X L lattice with periodic boundary conditions, i.e. (TL+i,y,z = 
ckx,L+i,z = (^x,i,z, <^x,y,L+i = (^x,y,i, for Convenience in the sequel. 
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A hint about the nature of the magnetic ordering in the model comes from 
considering the strongly anisotropic limit where = 0. In this case the horizontal 
plaquettes have zero coupling, which Hashizume and Suzuki gave the apt name of the 
“fuki-nuke” (“no-ceiling” in Japanese) model [HI[9]. The low-temperature order in such 
an anisotropic 3d plaquette Hamiltonian at = 0, 

L L L 


H, 


fuki-nuke — 


= - A 


x=l y=l z=l 


^x,y,z^x,y-\-l,z^Xyy+l,z-\-l^x,y,z-\-l 


L L L 




^x,y,z^x-\-l^y,z^x-\-l,y,z-\-l^x,y^z-\-l ? 


(4) 


x=l y=l z=l 

may be discerned by rewriting it as a stack of 2d nearest-neighbour Ising models. This 
can be carried out by defining bond spin variables Tx,y,z = o'x,y,zO'x,y,z+i at each vertical 
lattice bond. The r and a spins are thus related by 


cr. 


x,y,z ^x,y,l ^ D,y,2 '^x,y,3 ' ' ' D,y,z—1 i 


(5) 


where to maintain a one-to-one correspondence between the a and r spin configurations 
on a periodic lattice the constraints of the r-spins, flfli '^x,y,z = 1, must be satisfied 
and respected by an additional factor of 2 ^^^ in the partition function m- The resulting 
Hamiltonian when = 1 is then 


,) 


( 6 ) 


{'^x,y,zT'x+l,y,z + 'Tx,y,z'^x,y+l,z 

x=l y=l z=l 

which can be seen to be that of a stack of decoupled 2d Ising layers with nearest- 
neighbour in-layer interactions in the horizontal planes apart from the constraints that 
are expected to vanish in the thermodynamic limit m 

Since each 2d Ising layer will magnetize independently at the 2d Ising model 
transition temperature a suitable order parameter in a single layer is the standard Ising 
magnetization, which in terms of the original a spins is 

L L 


1 


T^2d,z — 


EE 


1^2 / . / ,^0!^,y,z^x:,y,z+l 
x=l y=l 


(7) 


The suggestion of Hashizume and Suzuki [8] and Johnston [12] was that similar order 
parameters could still be viable for the isotropic plaquette action. To avoid accidental 
cancellations between different planes we could use either absolute values for each plane 
or square the values. This gives a candidate 3d order parameter in the isotropic case as 
either 

L L L 

E EE-.. z'^x+l,y,z\ ) 5 (8) 

x=l y=l z=l 

or 


"^abs = 


1 

I? 


rrisq = 


L5 


E EE ^x,y,z^x-\-l,y,z 


(9) 


x=l \y=l z=l 




Macroscopic Degeneracy and order in the 3d plaquette Ising model 


4 


where we again assume periodic boundary conditions and with obvious similar 
definitions for the other directions, and For the isotropic model we 

would expect and similarly for the squared quantities, which can 

form a useful consistency test in simulations. 

Metropolis simulations indicated that and as defined above might 

indeed be suitable order parameters na for the isotropic plaquette model, but these 
were subject to the usual problems of simulating a strong first-order phase transition 
with such techniques. Here we discuss multicanonical Monte Carlo simulations [niE], 
combined with reweighting techniques [15], which allow us to carry out much higher 
precision measurements of and and confirm the suitability of the proposed 

order parameters. We also investigate their scaling properties near the first-order 
transition point. 

3. Numerical Investigation 

Details of the simulation techniques may be found in Ref. |6]. A two-step process 
is used, where estimates of an unknown weight function W{E) for configurations 
{cr} with system energy E = H{{a}) are iteratively improved. This replaces the 
Boltzmann weights that give the acceptance rate in traditional Metropolis Monte 
Carlo. In the first step the weights are adjusted so that the transition probabilities 
between configurations with different energies become constant, giving a fiat energy 
histogram [I6] as shown in Fig. below. The second step consists of a production 



e = E/L^ 


Figure 1. The approximately flat energy histogram p{E) for various lattice sizes 
produced in the multicanonical simulations on the left ordinate, and the reweighted 
canonical probability density at the temperature of equal peak-height on the right 
ordinate with dotted lines. 


run using the fixed weights produced iteratively in step one. This yields a time series 
of the energy, magnetization and the two different fuki-nuke observables ruabs and 
rUsq in the three possible different spatial orientations. With sufficient statistics such 
a time series together with the weights provides an 8-dimensional density of states 
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Q{E, m, ...) by counting occurrences of E, m, ... and weighting them with 
the inverse W~^{E) of the weights fixed prior to the production run. Practically, 
estimators of the microcanonical expectation values of observables are used. When 
measuring the order parameters, measurements are carried out every V = sweeps, 
because the lattice must be traversed once to measure the order parameters in all spatial 
orientations which has a considerable impact on simulation times. Skipping intermediate 
sweeps gives less statistics, but the resulting measurements are less correlated in the final 
time series. 

In Fig. 1^ we show the estimators of the microcanonical expectation values ((■)) of 
our magnetic observables O, 

{{o))(B) = ^on(£,o)/j;n(£,o), (lo) 

O O 

where the quantity VL{E^ O) is the number of states with energy E and value O of 
either the magnetization m or one of the fuki-nuke order parameters. We get an 
estimator for Vt{E, O) by simply counting the occurrences of the pairs {E, O) in the time 
series and weighting them with W~^{E). For clarity in the graphical representation in 
Figs. 1^ and we used a partition of 100 bins for the energy interval and an estimate 
for the statistical error of each bin was calculated by using Jackknife error analysis na 
with 20 blocks of the time series. 

The fuki-nuke parameters are capable of distinguishing ordered and disordered 
states, unlike the standard magnetization, which is revealed by Fig. It is also obvious 
from the first of Fig. that the different orientations of the fuki-nuke parameters 
are equal for the isotropic gonihedric Ising model, which we show for L = 20. This 
conhrms that the sampling is consistent in the simulation. We collect the microcanonical 
estimators for for several lattice sizes in the second of Fig. With the stored full 
time series and its weight function, we are able to measure the microcanonical estimators 
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e = E/L^ 


e = Eje 


Figure 2. Microcanonical estimators for the magnetization m and fuki-nuke 
parameter rUgq for lattices with size L = 20 and L — 27, where we used 100 bins 
for the energy normalized by system volume e = Ejl? for this representation. Errors 
are obtained from Jackknife error analysis with 20 blocks. 
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Figure 3. Left: Microcanonical estimators for the different orientations of the fuki- 
nuke parameters ^abs’sq ^ lattice with linear size L — 20, which fall onto two curves. 
The statistical errors are smaller than the data symbols and have been omitted for 
clarity. Right: Microcanonical estimators for the fuki-nuke parameter for several 
lattice sizes. 


for arbitrary functions of the measured observables /(O), 

((/(0)»(B) = Y, /(O) a(E,o)/Y»iE,o), (11) 

o o 

which provides a convenient way of calculating higher-order moments. For canonical 
simulations reweighting techniques |T5] allow system properties to be obtained in 
a narrow range around the simulation temperature whose width and accuracy are 
determined by the available statistics of typical configurations at the temperature of 
interest. Since multicanonical simulations yield histograms with statistics covering a 
broad range of energies (cf. Fig. it is possible to reweight to a correspondingly broad 
range of temperatures. The canonical estimator at hnite inverse temperature /3 > 0 is 
thus obtained as 

(0)(/3) = 5^{{0))(i5)e-y5^e-'’'^, (12) 

E E 

and Jackknife error analysis is again employed for an estimate of the statistical error. 

The behaviour of and in a Metropolis simulation [T2] is reproduced by 
the multicanonical data here as shown in Fig. Sharp jumps are found near the 
inverse transition temperature, as expected for an order parameter at a hrst-order 
phase transition. The peaks in the associated susceptibilities y = var(m) provide 
a suitable estimate of the finite-size pseudo-transition point. We find that the peak 
locations for the different lattice sizes L are fitted best by the modified first-order scaling 
laws with a leading 1/L^ correction appropriate for macroscopically degenerate systems 
discussed in detail in Refs. [S] and [B], 

/3^{L) = 0.551 37(3) - 2A6{3)/L‘^ + 2.4(3)/L^ , (13) 

where smaller lattices were systematically omitted until a fit with quality-of-fit 
parameter Q bigger than 0.5 was found. The fits presented have a goodness-of-fit 
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Figure 4. Canonical curves for the fuki-nuke parameters and over a broad 
range of inverse temperature /? for several lattice sizes L with the insets magnifying 
the sharp jump at the ordered phase. 



L-^ 

Figure 5. Inverse temperature /3{L) as determined from the peaks in y plotted against 
1/D (as appropriate for the modified scaling expected with macroscopic degeneracy). 


parameter Q = 0.64 and 12 degrees of freedom left. Fits to the other directions 
and hts to the peak location of the susceptibilities of give the same parameters 

within error bars and are of comparable quality. 

The inverse temperature /9(T) for the peak locations of the susceptibilities for 
both and is plotted against 1/L^ in Fig. [^together with the best £t curve 
giving the quoted values for the scaling coefficients. The estimate of the phase 
transition temperature obtained here from the finite-size scaling of the fuki-nuke 
order parameter(s), = 0.55137(3), is in good agreement with the earlier estimate 
= 0.551 334(8) reported in Ref. [6] using fits to the peak location of Binder’s energy 
cumulant and specific heat and the value of (3 where the energy probability density has 
two peaks of the same height or weight. 

As a visual confirmation of the scaling, the normalized susceptibilities for 

various lattice sizes are plotted against f3 + 2.46/L‘^ — 2.4/D (where the shifts are 
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Figure 6. Canonical curves for the susceptibilities associated with the fuki-nuke 
parameters and plotted against the shifted inverse temperature /3+2A6/D — 
2A/D for several lattice sizes L. The peaks fall onto the same point. 


determined by the fitted scaling corrections) in Fig. It is clear that the peak positions 
fall on the same point. In Fig. we plot the peak values of both Xmi^^ and Xmi^ divided 
by the system volume against 1/L. Empirically the plotted hts for the maximum values 
are 

{L)/L^ = 0.135 20(29) - 0.757(10)/L + 1.33{7)/L^ (14) 

^abs 

where we used all lattices with sizes greater than or equal to Lmin = 10 (with 15 degrees 



Figure 7. Maximum value of the normalized susceptibility plotted against 1/L 

for both and . 
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of freedom left and a goodness-of-fit parameter of Q = 0.07) and 

= 0.134 86(14) - 0.174(4)/L - 1.413(27)/L2 (15) 

sq 

with Lmin = 8 (with 17 degrees of freedom and Q = 0.08). The leading corrections of 
the peak valne of the specific heat divided by the system volume are of order 0(1/L^), 
so those for the susceptibility are already much stronger. Forcing a fit to leading 1/L^ 
corrections here gives slightly poorer fits and we cannot distinguish empirically. The 
constants, however, are stable around 0.135, and barely change with the 

various fits, showing that the proportionality factor of oc is independent of the 
leading corrections. By analogy with the maximum of the specific heat, 13 EllI] we 
expect for the extremal values of the susceptibilities 

= r j +..., ( 16 ) 

where the gaps Afriabs.sq = '^abs)sq'^ ~ of fuki-nuke order parameters in 

the inhnite system enter. This is consistent with the measurements since a value of 
^ 0.135 implies Arriabs.sq ~ 0.99, which is plausible from the insets in Fig. g 
giving an impression of the value of as L —)■ cx) and noting that = 0. 

4. Conclusions 

The macroscopic degeneracy of the low-temperature phase in the 3d plaquette gonihedric 
Ising model excludes standard magnetic ordering. However, consideration of the strongly 
anisotropic limit of the model suggests that a planar, fuki-nuke order may still be 
present. Multicanonical simulations of the model strongly support this suggestion, with 
the various fuki-nuke magnetizations all showing order parameter like behaviour. If 
the effects of the macroscopic low-temperature phase degeneracy on the corrections to 
scaling detailed in Refs. [5] and [6] are taken into account, the estimates for the transition 
point obtained from these fuki-nuke magnetizations are fully consistent with estimates 
obtained from energetic quantities. There are, however, stronger finite-size corrections 
to the peak value of the susceptibilities for the fuki-nuke order parameters. 

Acknowledgements 

This work was supported by the Deutsche Forschungsgemeinschaft (DFG) through 
the Collaborative Research Centre SFB/TRR 102 (project B04) and by the Deutsch- 
Franzosische Hochschule (DFH-UFA) under Crant No. CDFA-02-07. 

[1] G. K. Savvidy and F. J. Wegner, Nucl. Phys. B 413 (1994) 605; 

G. K. Savvidy and K. G. Savvidy, Phys. Lett. B 324 (1994) 72; 

G. K. Savvidy and K. G. Savvidy, Phys. Lett. B 337 (1994) 333; 

G. K. Bathas, E. Floratos, G. K. Savvidy, and K. G. Savvidy, Mod. Phys. Lett. A 10 (1995) 
2695; 

G. K. Savvidy, K. G. Savvidy, and F. J. Wegner, Nucl. Phys. B 443 (1995) 565; 

G. K. Savvidy and K. G. Savvidy, Mod. Phys. Lett. A 11 (1996) 1379; 



Macroscopic Degeneracy and order in the 3d plaquette Ising model 10 

G. K. Savvidy, K. G. Savvidy, and P. G. Savvidy, Phys. Lett. A 221 (1996) 233; 

D. Johnston and R. P. K. C. Malmini, Phys. Lett. B 378 (1996) 87; 

G. Koutsoumbas, G. K. Savvidy, and K. G. Savvidy, Phys. Lett. B 410 (1997) 241; 

J. Ambj0rn, G. Koutsoumbas, and G. K. Savvidy, Europhys. Lett. 46 (1999) 319; 

G. Koutsoumbas and G. K. Savvidy, Mod. Phys. Lett. A 17 (2002) 751. 

[2] A. Lipowski, J. Phys. A: Math. Gen. 30 (1997) 7365. 

[3] D. A. Johnston, A. Lipowski, and R. P. K. G. Malmini, in Rugged Free Energy Landscapes: 

Common Computational Approaches to Spin Glasses, Structural Glasses and Biological 
Macromolecules, ed. W. Janke, Lecture Notes in Physics 736 (Springer, Berlin, 2008), p. 173. 

[4] R. Pietig and F. Wegner, Nucl. Phys. B 466 (1996) 513; 

R. Pietig and F. Wegner, Nucl. Phys. B 525 (1998) 549. 

[5] M. Mueller, W. Janke, and D. A. Johnston, Phys. Rev. Lett. 112 (2014) 200601. 

[6] M. Mueller, D. A. Johnston, and W. Janke, Nucl. Phys. B 888 (2014) 214. 

[7] M. Mueller, W. Janke, and D. A. Johnston, Physics Procedia 57 (2014) 68. 

[8] Y. Hashizume and M. Suzuki, Int. J. Mod. Phys. B 25 (2011) 73; 

Y. Hashizume and M. Suzuki, Int. J. Mod. Phys. B 25 (2011) 3529. 

[9] M. Suzuki, Phys. Rev. Lett. 28 (1972) 507. 

[10] M. Mueller, W. Janke, and D. A. Johnston, Nucl. Phys. B 894 (2015) 1. 

[11] C. Castelnovo, C. Ghamon, and D. Sherrington, Phys. Rev. B 81 (2010) 184303. 

[12] D. A. Johnston, J. Phys. A: Math. Theor. 45, 405001 (2012). 

[13] B. A. Berg and T. Neuhaus, Phys. Lett. B 267 (1991) 249; 

B. A. Berg and T. Neuhaus, Phys. Rev. Lett. 68 (1992) 9. 

[14] W. Janke, Int. J. Mod. Phys. (7 03 (1992) 1137. 

[15] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 61 (1988) 2635; 

W. Janke, in Computational Physics: Selected Methods - Simple Exercises - Serious 
Applications, eds. K. H. Hoffmann and M. Schreiber (Springer, Berlin, 1996), p. 10. 

[16] W. Janke, Physica A 254 (1998) 164; 

W. Janke, in Computer Simulations of Surfaces and Interfaces, eds. B. Diinweg, D. P. Landau, 
and A. 1. Milchev, NATO Science Series, H. Math. Phys. Ghem. 114 (Kluwer, Dordrecht, 2003), 
p. 137. 

[17] B. Efron, The Jackknife, the Bootstrap and other Resampling Plans (Society for Industrial and 

Applied Mathematics, Philadelphia, 1982); 

W. Janke, in Proceedings of the Euro Winter School Quantum Simulations of Complex Many- 
Body Systems: From Theory to Algorithms, eds. J. Grotendorst, D. Marx, and A. Muramatsu, 
NIG Series Vol. 10 (John von Neumann Institute for Computing, Jiilich, 2002), p. 423; 

M. Weigel and W. Janke, Phys. Rev. Lett. 102 (2009) 100601; 

M. Weigel and W. Janke, Phys. Rev. E 81 (2010) 066701. 



